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Set Books 
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It is essential to have these books; the course is based on them and will 


not make sense without them. 


Conventions 


Before working through this correspondence text make sure you have read 
A Guide to the Linear Mathematics Course. Of the typographical conven- 
tions given in the Guide the following are the most important. 


The set books are referred to as: 


K for An Introduction to Linear Analysis 
N for Linear Algebra and Matrix Theory 


All starred items in the summaries are examinable. 


References to the Open University Mathematics Foundation Course Units 
(The Open University Press, 1971) take the form Unit M/00 3, Operations 
and Morphisms. ` 
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9.0 INTRODUCTION 


In this unit we return to the study of linear differential equations, which we 
started in Unit 4, Differential Equations I. There we looked at linear dif- 
ferential equations from a general point of view, and saw how to reduce a 
very general type of first-order differential equation to the evaluation of 
integrals. For equations of higher order, there is no such general method 
for obtaining analytic solutions, but there are methods for dealing with 
special types of equation which often arise in applied mathematics. There 
are also general results yielding information which can be very useful, 
even though it may fall short of a complete solution. 


The most fundamental of these results is the one we have seen already; 
that is, that the problem of finding the general solution of a linear differen- 
tial equation of the form 


Ly=h 
where h is a known function and L is a linear differential operator, can be 


split into two parts. One part is to find the general solution y, of the asso- 
ciated homogeneous equation 


Ly =0; 
the second is to find a particular solution y, of the original equation. Then 


it follows from the general theory of linear problems that y, + y, is the 
general solution of the original equation. 


In the present unit, we tackle the first of these parts: finding the general 
solution of a homogeneous differential equation of the form Ly =0, 
where Z is a linear differential operator of order at least 2. We shall deal 
first with a special, but very useful, type of homogeneous equation for 
which there is a recipe for writing down the general solution in all cases. 
These are the equations with constant coefficients. Afterwards we shall deal 
with homogeneous equations in general. Nonhomogeneous equations will 
be dealt with in Unit 11, Differential Equations III: Nonhomogeneous 
Equations. 


The key idea is the fact that the solution set of Ly = 0 is a vector space 
whose dimension is equal to the order of L. This has the consequence that 
if we can find n linearly independent solutions of Ly = 0 (where n is the 
order of L) then these » solutions form a basis for the solution space. The 
entire solution space is therefore the set spanned by these n solutions, and 
the general solution is an arbitrary linear combination of them. We saw an 
example of this in the Foundation Course (Unit M100 31, Differential 
Equations II), where we studied the differential equation 


X'()- X(n-0 — (te[0, oo) 


which describes the free vibrations of a mass-spring system with suitable 
values for the mass and the stiffness of the spring. 


X(t) 
In the notation of this course (using y in place of X), the equation can be 
written 
y'*y-0 
or 
(D? + ly 20 
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where D is the differential operator. Since the differential operator D? + I 
is of second order, the solution space has two dimensions. Two linearly 
independent solutions of the equation are the sine and cosine functions, 
and since they are linearly independent, they form a basis for the solution 
space. The solution space is therefore 


(cos, sin» 
and the general solution (an arbitrary element of this space) is the function 


€, COS +c, Sin 


t~ ¢, cost+e,sint (te [0, œ), 
where c, and c; are arbitrary numbers, the so-called arbitrary constants. 


Whatever the differential equation, provided it is linear and homogeneous, 
the method of approach is the same: find » different solutions, check 
that they are linearly independent, and construct the general solution from 
them. If the equation has constant coefficients (like the equation y" + y 0 
mentioned above) then it is always possible to find the necessary n solutions 
and hence solve the equation: this is discussed in Section 1 of this unit. If 
the equation does not have constant coefficients (as, for example, in the 
equation y" + xy — 0), and the order is 2 or greater, then it is not, in 
general, possible to find explicit solutions; but it is still possible to find out 
a lot about the solutions. In Section 2 of the unit we Shall see, in parti- 
cular, how to prove the principle theoretical result for these equations, 
that the dimension of the solution space is equal to the order of the equa- 
tion; we also obtain a method for checking that a given set of n solutions 
of the equation really are linearly independent and therefore can be used as 
a basis. 


The first section of the unit is based on Sections 4-1, 4-2 and 4-3 of K; the 
second is based on Sections 3-4, 3-5 and 3-6. We reverse K’s order of 
presentation in order to put the easier material first, but if you prefer to 
treat general theorems before particular cases, there is no objection to 
your studying the second section of the unit before the first. 
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9.4 CONSTANT-COEFFICIENT EQUATIONS 
9.1.1 Polynomials in D 


As we saw in Unit 4, Differential Equations I, the general homogeneous 
linear differential equation has the form 


a oy" o) + a, -,Q)y 7 Qe + 


*aC)yQ)-ra()y0)-0 (ye CN) 
where J is some interval of R and a,, a,.,, ..., a, ag are continuous real 
functions with domain J, and a, is not the zero function. The functions 
an» 05-1, +++, dg are called the coefficients, and if they are constant func- 
tions, then the equation is said to have constant coefficients. The signifi- 
cance of the fact that the coefficients are constants is indicated in the first 
reading passage. 


READ page K126 as far as “. . . with relative ease". 


Notes 


(i) Equation (4-1) In this unit we shall be concerned only with the case h — 0. 
(i) line 4 Only a, has to be nonzero (if it were 0, the order of the equation 
would be less than 7). 
(ii) lie 1] “normal form" simply means that we divide the whole equation 
by a, (which we know is nonzero) and then re-define an-ı, ..., ao to stand 
for the quantities previously represented by an-1/ān, ..., @o/@n- 
(iv) line —9 “Algebraically such operators . . ." The manipulation of these 
` operators is considered in more detail in Example 4 on pages K51-52. 


The relationship between constant-coefficient linear operators and poly- 
nomials is very important so we will go into it in a little more detail than K. 
Given such a linear operator 


L= D" + ap- D"! + + ao 
we associate with it a polynomial P such that 
P(x) =x +a H Hao (xER) 


In this way we obtain a mapping M from the set of constant-coefficient 
linear differential operators (CCLO) with addition, scalar multiplication 
and the (composition) product to the set of polynomials with formal sums, 
scalar products and products, which is an isomorphism. 


M 
CCLO 4——2 Polynomials 
M^! 


For the product of differential operators in general we cannot “ multiply 
out” as we do for polynomials, as can be seen in the example on page K88: 


(xD + 2)(2xD + 1) = 2x? D? + 7xD +2 
whereas straightforward formal polynomial multiplication wouid give 
2x? D? + 5xD +2 


If Lı, Lz are constant-coefficient operators the situation is much nicer. If 
the corresponding polynomials are P,, P}, we can write 


L, = M(P), L; = M7*(P2) 


and then L, + L, corresponds to the polynomial P, + P2, ie. Ly +L, = 
M-P, + P3); 


the operator aL, corresponds to the polynomial aP,, i.e. aL, = M~'(aP,); 


and the operator L; ° L =L,L corresponds to P,P}, ie. LL, = 
M~\(P,P,). 
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Hence we have the following commutative diagrams (indicating the 
presence of a morphism) 


M ; 
CCLO— — — — —9 Polynomials 


scalar 
multiplication 
of operators 
by reals 


scalar multiplication 
xa| of polynomials 
by reals 


M " 
CCLO— — — — —» Polynomials 


and the two diagrams 


M 
CCLO x CCLO— — — — —^ Polynomials x Polynomials 
addition of addition of 
+ + b 
operators polynomials 


+ M Y 
CCLO————————>Polynomials 


M 
CCLO x CCLO— — — — —» Polynomials x Polynomials 


composition 
product of 
operators 


product of 
polynomials 


Y M 
CCLO— — — —— —» Polynomials 


For example 


D? + 1—— —— —93* 41 


x2 x2 


M 


2D? 4:2 — — oaa 


D +1,D-1 xd 2 
> SE LL 


+ + 
M 
D* + D—— —— —À +x 
since 
KD? + 1) + (0 = D]y = (D? + Dy + (D = Dy 
= D'y y Dy-y 
= D?y + Dy 
=(D? + Diy 
Exercise 


If L, and L, are constant-coefficient operators, show that L, o L, = L, o L 
by appealing to the mapping : 


M: CCLO ——- Polynomials 
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Solution 


Since M is an isomorphism of (CCLO, ») to (Polynomials, x ) and 
multiplication of polynomials is commutative, it follows that 
L, ° L; =L,°L,. 


Exercise 
Draw the commutative diagram associated with the (composition) product 
of (D? + 1) and (D — 1). 


Solution 
D? +1, D— 1———>x? + 1, x-1 
o x 
D-D.p-l = x) x 4x1 
Since 


I + DO — Diy = (P? + DEC = y] 
= (D? + 1)(Dy — y) 
= D'(Dy — y) + Dy - y 
= Dy- D?y+ Dy -y 
=(D? — D? + D — l)y 
For the reasons mentioned in the reading passage, it is important in this 
kind of work to know how to factorize a polynomial into linear and 
quadratic factors. 
Example 
Factorize the operators D? — 8, D? — 8. 


As explained in the text following note (iv), the problem is the same as the 
problem of factorizing the polynomials x? — 8, x? — 8. Now any poly- 
nomial P(x) of degree n can be factorized into the form 


P(X) = a(x — e4)(x — 85) +>: (x — 9) (1) 


where a, is the coefficient of x" and a4, 2;,..., «, are real or complex 
numbers called the roots of P(x) — 0. They can be found by using the fact 
that every root is a solution of P(x) — 0. 


For the polynomial x? — 8, we solve 
x? -8=0. 
The roots are + J8, so the factorization of x? — 8 is 
x? — 8 = (x — 8x + 8) 
and the factorization of D? — 8 is (D — JID + /8). 
For x? — 8, we must solve 
x -820. 


One root is obviously x — 3/8 — 2, and so there is a factor x — 2. To find 
the others we try 


x? — 8 = (x — 2)(ax? + bx +) 
(there are two more factors; so their product will be a quadratic). 
Multiplying out the right-hand side gives 


x? — 8 = ax? + (b — 2a)x? + (c — 2b)x — 2c 


We can regard the x?, x?, x and 1 as a basis for the vector space P,, of 
polynomials of degree three or less. Then the two sides of the equation are 
apparently two different ways of representing the same polynomial in 
terms of this basis. But we know that any element of a vector space can be 
expressed uniquely in terms of a basis: so we can equate corresponding 
coordinates (“equate coefficients of like powers of x”, as it is usually put in 
this context), i.e. 


l=a 
0=b-2a 
0=c -2b 
—-8 = -—2c 


Solving these in succession gives a — 1, b — 2, c — 4 and again c — 4, and 
putting these numbers into ax? + bx +c in the previous factorization 
gives 


x? — 8 = (x —2)(x? + 2x + 4) (2) 
To complete the factorization, we would solve 
x? 42x+4=0. 


If we apply the formula for quadratics to this, however, we find that both 
the roots are complex; thus the linear factors of x? +2x+4 are also 
complex. Such factors do not interest us here, however; they would not 
lead to a suitable factorization of D? — 8, because our definition of a 
linear differential operator (see page K126, line 4) requires the coefficients 
to be real. Thus, the simplest factors we can get for D? — 8 are the ones 
corresponding to Equation (2), namely 


D? — 8 = (D —2YD* + 2D +4). 


Exercises 


1. Exercises 13(a), (b), (c), (d), page K90. In (d) one root of the relevant 
polynomial equation is 2. *Irreducible" means that any further 
factorization would involve complex coefficients. 

Exercises 5(a), (b), (d), page K127. In each case, at least one of the 
roots of the corresponding polynomial equation is an integer. Use 
the fact that any integer root of 


S 


X RaQ2aQx'"i-ecRaxaay-0 


must be a (positive or negative) factor of ag. 


Solutions 


1. (a) D'—-3D42-(D—1Y(D-2) 
(b 2D*-5D--2-2X(D 4 3D -2) or 2D + 1)(D + 2) 
(c) 4D? +4D4+1=4(D +4) or (2D + 1)? 
(d D -3D +4=(D- 2D +1) 


One root of x? — 3x? + 4 = 0 is 2, so one factor is x — 2. 


Trying x? — 3x? + 4 = (x — 2Y(ax? + bx + ce) 
= ax! + (b —2a)x? + (c — 2b)x — 2c 


and equating coefficients of like powers of x gives 


g—l;b—2a- 3;¢-2b=0; —2c — 4. 
Thus 
a=1,b=-le=-2,¢= —2 


again. 
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Hence, the other factor is (x? — x — 2), which factorizes to 
(x — 2)(x + 1) and we get 


x — 3x7 + 4 = (x — 2)(x — 2)(x + 1). 
2. (a) D'+4D?+5D4+2=(D+ 1)(D+2) 


One root of x? + 4x? + 5x - 2-0 is —1, so x+ l isa 
factor. 


By the method used in Solution 1(d), the other factor 
is x? + 3x + 2, which factorizes to (x + 1)(x + 2) and we 
get 


x 44x? + 5x +2 = (x + L(x + 1x + 2) 
(b D?—D? 4 D—1=(D— 1D? 4 1) 


The quadratic x? + 1 = 0 has no real roots, so there are 
no simpler real factors. 


(d D*—5D!4-(D-— YD 4 D(D —2Y0D + 2). 
One root of xt — 5x? + 4 — 0 is 1, so a factor is x — 1. 
This gives 
xt — 5x? + 4 = (x — DG +x? — Ax — 4). 
One root of x? + x? — 4x —4=0 is —1, and so 


X +x? — 4x — 4 = (x + Ix? — 4) 
=(x + 1) — 2 (x + 2) 


9.1.2 Second-order Equations 


The next reading passage explains how the factorization technique of the 
preceding section enables us to solve linear constant-coefficient equations. 


READ the whole of Section 4-2, starting on page K127. 


Notes 


(i) line —3, page K127 A lemma is a preliminary minor theorem. “ Null 
space" means “kernel”. 

(i) line 1, page K128 In other words, we prove that L,y —0 implies 
(Lı ...L)y =0, for i—1,...,n. As K remarks the proof is trivial: 
almost all that is being used is the commutativity of operator multiplica- 
tion to get the L; to the end. 


solution set of solution set of 
Liy=0 Lı Lo. ..L,. y-0 
solution set of 
L2y=0 


(iii) 


(iv) 


(v) 


(vi) 


(vii) 


(viii) 
(ix) 


(x) 


line 12, page K128 The equation (D — 2)y =0, i.e. y’ — 2y — 0, could be 
solved by the method of separation of variables described on page K96; 
but it is enough just to check by differentiation that y: x ———9 e?* is a 
solution. 

line 13, page K128 The Wronskian is a method of testing sets of func- 
tions for linear independence, which we shall study later in the unit. 
An alternative method that works here is this: if cie?*+c2e7?* is 
the zero function, then taking the particular values x =0 and x =1, we 
have c; + c; — 0 and c,e? + e; e7? = 0, and it follows (since e? # e^?) 
that e; = c; =0. 

lines 3 to 5, page K129 Do not concern yourself here with the method 
by which the second solution y(x) = xe™ is obtained; a method for obtain- 
ing this solution is given in Unit 11, Differential Equations III. The 
important point here is that if the auxiliary equation has equal roots, then 
a second solution of the differential equation is x times the first. 

line 9, page K129 The reason why the method breaks down is that the 
e*!* and e??? are now complex; consequently the functions 


x > ec *bx 
and 
x > e- bi)x 


are not real functions, because their codomains are sets of complex 
numbers, rather than 2. Consequently, these functions do not belong to 
the domain of the linear operator D? + aD + a; associated with the 
equation we are considering; for this domain is defined (on page K86) to 
be (i), which is a set of real-valued functions, (i.e. functions with 
codomain 4»). Of course we could have defined this domain to include 
complex-valued as well as real-valued functions, and had we done so 
functions such as e'^*^^* would have been perfectly acceptable solutions 
of the equation. (There would still remain the problem, of course, of what 
we mean by differentiating such a function: but this can also be overcome 
reasonably easily.) It is often convenient, in fact, to extend the domain to 
include complex-valued functions and so take advantage of the greater 
simplicity of complex exponentials such as e+’) in comparison with 
expressions such as e° cos bx. For the present, however, we stick to real 
functions in our rigorous arguments, and treat the argument based on 
e+ as a piece of heuristic reasoning (i.e. plausible reasoning) which 
must be checked afterwards by rigorous methods. (There is a discussion, 
similar to that of Case 3 in K, in Unit M100 31, Differential Equations II.) 
line 16, page K129  Euler's formula was discussed in Unit M100 29, 
Complex Numbers II. Note that the symbol “x” in this formula corres- 
ponds to the “bx” of the rest of the calculation. 

line —12, page K129 Here cs and c4 are defined as c; + c; and i(c, — c2), 
respectively. 

line —11, page K129 Mathematicians often use the word “formal” to 
mean that symbols are being manipulated without adequate regard for 
their meaning, i.e. for the rules governing these manipulations. 

last line of table, page K130 There is a small change of notation: the 
numbers denoted by c; and cs on the preceding page are denoted here by 
c, and c2. 


Examples 


Find the general solution of the following differential equations 
@ (D'-2D-3y-0. 

(b (D? +4D + 4)y=0 

(c) (D? + 6D + 10)y =0 


Solutions 


12 


(a) We see that the operator factorizes 
(D? -2D—3)- (D - 3(D +1) 


and hence the auxiliary equation has distinct real roots 3 
and —1. This gives us the general solution 


W(x) = e&* + 0, e7*, 
(b) In this case the operator factorizes 
(D? + 4D + 4) =(D + 2) 
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(c) 


Exercises 


l. Exercis 
tution i 
2. Exercis 
Unit 7, 


Solutions 
I. 


to give the repeated root —2 of the auxiliary equation. 
Hence, we have as a general solution 


WR) = (e + c3 x)e7?*-. 
In this case we cannot find real factors of the operator. We 
can use the formula to calculate the roots of the auxiliary 
polynomial equation: 
—6 +,/36 — 40 
2 
or we can observe that the operator is of the form D? — 
2aD +a? + b? where a = —3 and b= 1. In either way we 


see that the roots of the auxiliary equation are —3 +i and 
the general solution is 


=-34i 


Wx) = e4e7?* cos x + c; e^? sin x. 
The fact that an operator of the form D? — 2aD + a? + b? 


corresponds to solutions of the form e™ cos bx, e^" sin bx 
will be useful to remember for later exercises. 


es 1, 4, 9 and 14, page K130. Check your solutions by substi- 
n the differential equation. 

es 18 and 20, page K130. (We met initial-value problems in 
Recurrence Relations.) 


Kl. Differential equation. (D^ + D—2)y «0 


(DE) 
Auxiliary equation — m? +m—2=0 
(AE) 
Solutions of AE m=1 or -2 
Solution of DE Wx) = ce" + eze~? 


K4. Differential equation (D? —2D)y =0 

Auxiliary equation — m? — 2m =0 

Solutions of AE m-0or2 

Solution of DE y(x) = cy + c; e?* (since e? = 1) 
K9. Differential equation. (D? - 2D + 2)y - 0 

Auxiliary equation — m? — 2m - 220 

Solutions of AE m-lctil-i 

Ge acl b = 1) 
Solution of DE y(x) = e*(c, cos x + c; sin x) 


Ki4. Differential equation (9D? -- 6D -- 1)y «0 


Auxiliary equation — 9m? + 6n 4-120 


Solutions of AE m = —1 (twice) 
Solution of DE Wx) = (e, + e; x) exp (33). 
2. K18. Differential equation (4D? — 12D --9)y 20 

Solutions of AE m = } (twice) 
General solution 

of DE is y = (e, + e; x) exp (3x) 
Conditions y0) = ¢, = 1 

yO =e, +a =F 

whence € =1,¢,=2 
Required solution 

of DE y(x) = (1 + 2x) exp (3x) 

K20. Differential equation (4D? — 4D + 5y =0 

Solutions of AE m=tti 
General solution 

of DE y = eU? (c cos x + c, sin x) 


13 


Conditions 30) =e, =4 
YO 2 $e +e! 
whence o-Lo-i 
Required solution 
of DE y(x) = et?» cos x + $ sin x) 


9.1.3 Equations of Order Greater than 2 


The method we have been using for second-order homogeneous equations 
Beneralizes in a straightforward way to equations of higher order. The 
method is described in the next reading passage. 


READ the whole of Section 4-3, starting on page K132 omitting the proof 
of Theorem 4-1. 


Notes 


(i) line 15, page K132 “these functions are linearly independent" Take this 
for granted at present; the authors return to the question of linear inde- 
pendence in Example 5 on page K135. 

(i) line —8, page K135 As we have already remarked, the Wronskian is a 
method for testing linear independence of functions, which we shall look 
at later in the unit (sub-section 9.2.4). 

(iii) Equation (4-14), page K135 This is the usual strategy for proving linear 
independence: weshow that the only linear combination of these functions 
that is equal to the zero function is the combination with €, 0m 
£s =0. 

(iv) line —2, page K135 Thesymbol = (read “ identically equals”) is a way of 
indicating that this is a relation between functions, not numbers: the func- 
tion ¢:(D + 5))(D? — 4D + 13)e?* is equal to the zero function. K does 
not use this consistently: for instance, it could also be used in Equation 
(4-14). 

(v) line 13, page K136 To prove that cs and cs are 0, use the fact that e?* is 
never 0, and consider the particular values x = 0 and x = iv in the equation 
cs COS 3x + ce sin 3x — 0. ` 


Exercises 


l. Exercises 1,3,11 and 14, page K136. 

2. Exercises 17(a) to (d), page K136. (Hint The necessary formulas are 
at the top of page K134.) 

3. Find the solution of y" + 3y” — y — 3y — 0 that satisfies the initial 
conditions y(0) = 2, y'(0) = —2, y"(0) = 10. 

Solutions 


1. Kl. The DE is (D + 3D? — D — 3)y — 0. Factorization by 
the method of sub-section 9.1.1 converts this to 


(D +3XD + 1(D— ly =0 
The general solution is therefore 


YX) = ae? 


"toe oe. 
K3. The DE is (4D? + 12D? + 9D)y — 0, ie. 
4D(D + 33 y — 0. 
The general solution is therefore 
W(X) = e + (c2 + c3 x) exp (—4). 
Kil. The DE is (D* + D? + Dy =0, i.e. 
DD? + D 4 1y —0. 


The roots of m? + m + 1 =0 are ERIS 


14 
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The general solution of the DE is 
JQ) =e, + cx 
+e UP (c, cos 4 NEUES ca sin 3 /3x) 
K14. The DE is (D5 + 2D? + D)y — 0, which factorizes to 
D(D* + 2D? + y =0 
or D(D? + 1)?y =0. 
The general solution is therefore 
(x) = ey + (e; + c3 x) cos x 
+ (c4 + c5 x) sin x. 
2. KI7. (a) Since x?e*"! = exe", a suitable operator is 
(D — 1} or (D — 1} times any polynomial in D. 

(b) A suitable operator is D? — 2a D + (a? + b?) with 
à —2, b 2, ie. D* — 4D + 8. (See example on 
page C12.) 

(c) (D? +1)? annihilates both 2x? sin x and x sin x, 
and hence also their sum (or in fact any function 
x ——> p(x)sin x, where p is a polynomial of 
degree 2 or less). 

(d) D? annihilates 3-- 4x and D 2 annihilates 
—2e7?*, Their product, D?(D +2) annihilates 
both and is a solution to the problem. 

3. By Part ] of Solution 1, the general solution of 


y"-3y'- y —-3y 20 


y(x) = ee7 9 + ee roe 
The conditions are y(0) — c, +c +c, =2 
y(0)2 —3¢, — c; + c3 = —2 
y'(Q =9e, +e, +e, —10 


The solution of these equations (use Gauss elimination) is 
€, =1, c5 =0, c4 — 1, and so the required solution of the 
differential equation is y(x) = e^?* + e. 


Damping and Stability of Vibrations 


Linear constant-coefficient equations of the second order are very useful 
in the theory of mechanical and electrical vibrations. We saw in Unit 
M100 31, Differential Equations II, how the differential equation 


mX'(t) + sX(t) 20 


with m and s positive numbers, provided a model of a mass-spring system. 
This equation can be made more realistic by including a term to describe 
the frictional forces. If we assume that the frictional force is proportional 
to the velocity X'(1) and is in the opposite direction, then we can allow for 
this in the equation of motion by adding to the spring force —5X(1), an 
additional frictional force — k X'(t), where k is a positive number. 


motion 
OO 


frictional force OSH 
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Thus the equation of motion becomes 
mX"(t) = —sX(t) — kX'(t), 


which is a second-order linear constant-coefficient equation. It can be 
written in the equivalent form 


(nD? + kD + $)X =0. 


The auxiliary equation (using p for the polynomial variable instead of m, 
which is used for mass here) is 


mp? 4 kp s-0 
and this has roots which are 
—k X Jk? — 4ms 
2m ` 


The methods given in sub-section 9.1.2, and summarized on page K130, 
show that the solutions of the differential equation can be of three types: 


G) small damping: 0 < k? < 4ms. The auxiliary equation has complex 
roots. In this case the general solution is 


X(t) = exp(—kt/2m) 


JAms — k? 
2m 


t+ cy, sin 


( AJ Ams — à ) 
x {c cos —— —— (t 
2m 

Typical motions are indicated in the diagrams: the motion is oscil- 
latory, as in the undamped case k — 0, but now the amplitude de- 
creases with time. Mathematically, the oscillatory motion is provided 
by the cosine and sine terms, and the damping by the exponential 
term. 


X(t) X(t) 


(ii) critical damping: k? = 4ms. The auxiliary equation has equal roots, 
and the general solution is 
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x(t) + X(t) 


> 


t 


This time the motion is not oscillatory. 
(iii) large damping: k? > 4ms. The auxiliary equation has real roots, 
and the general solution is 


—kt 
X(t) = exp (=) 
Jk? — 4m 
2m 


s 
t+ c, exp 


x (s exp 
Again the motion is not oscillatory. 


4X() ^x(t) 


» 


t 


In each case the character of the motion depends crucially on the fact that 
k is positive. If k were negative, the factor exp (—kt/2m) would increase 
instead of decrease, and the motion would get more and more violent in- 
stead of dying away. This may seem unrealistic, but it can happen, as the 
following (optional) example shows. 


Example (Optional)* 


Consider the stability of the vibrations of an aircraft wing. In the analysis, 
we shall assume for simplicity that each section of the wing is elastically 
coupled, as if by a spring, to a fixed point (rather than to other parts of 
the wing). We also assume that the air moves past the wing in a horizontal 
direction (as in a wind tunnel) rather than the wing past the air. We denote 
the velocity of the air by V, and take it to be constant and positive. 


There are two vertical forces on the wing. One is the elastic force, which 
depends on X(t), the vertical displacement of the wing section from its 
equilibrium position. We denote this force by F(X(t)), with the convention 
that an upward force is positive, The second is the aerodynamic lift force, 
which depends on the angle of incidence a that the wing makes with the 
direction of the air flowing relative to it; we denote this lift force by L(a). 


* This example ends on page 20. 


a 
ri 


actual 
position 


equilibrium 
position 


Now, if the wing is moving upwards with velocity X’(), the motion affects 
the angle of incidence; for during the time a particle of air travels a distance 
V to the right, the aerofoil travels a distance X'(t) upwards and so the 
particle's displacement relative to the aerofoil is a combination of V to the 
right with X'(t) downwards. 


path of air particle x) 


relative to aerofoil 


Thus, the velocity of the air stream relative to the aerofoil makes an angle 0 
with the horizontal, where 


XO 
tan 0 = Ty 


For small 0 we can approximate* tan 0 by 0, so we have 0 ~ X'(t)/V, and 


the angle of incidence « is given (see diagram) by 


i 
su- 20 


where a is the angle between the aerofoil and the horizontal. 


* This is the Taylor a 


pproximation of degree | at0: see Unit M100 14, Sequences and 
Limits II. 
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Combining the two vertical forces on the aerofoil, we obtain the equation 
of motion 


mX"() = F(X()) + 1% = zO) f 


where m is the mass of the aerofoil. 


Let us make the assumption that F and L can be approximated by linear 
functions (i.e. polynomial functions of degree 1). We can use the Taylor 
approximation of degree 1 (see Unit M100 14, Sequences and Limits II), 
which for a general function f is 


S) = f (%0) + (x — xo) f). 
In our case x, is taken as the initial equilibrium position, when we have 
X(0) = X'(0) — 0. 


We then have 
F(X) = F(0) + X()F(0) 
and 
X'(t) 
y 


L(a - ZË) = neo - reo 


Substituting these results in the differential equation, we get 
X'() 


mX"(r) = F(0) + F'(0)X(t) + Lo) — L'(ao) Y 


Since we are taking X(f) to be the displacement from equilibrium, this equa- 
tion must have the zero function as a solution; that is, we must have 


0 = F(0) + L(X9) 
The equation of motion thus simplifies to 


X'(t) 
4 


mX"(t) = F(0)X(r) — L(G) 


which we shall write as 


mX"() + £o) X(t) + sX(t) 20 


where s = — F'(0); this quantity is positive since an upward displacement 
of the aerofoil leads to a downward force on it from the spring. 
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The simplified equation of motion we have just obtained is of the linear 
type we have already investigated. As we have seen in this sub-section, 


kt ; : " 
thesolutions contain a factor exp (- =) where k is the coefficient of X'(t) 
nm, 


in the differential equation. If this coefficient is positive, any oscillations 
will die away, but if it is negative, they will increase. Thus, the condition 
for the wing to be stable under this type of vibration is that k should be 
positive, i.e. (since V is positive) 


L(a) > 0 


The lift must increase, not decrease, with the angle of incidence. The func- 
tion L normally looks like that in the diagram. 


a 
L@) 


Y 


Thus stability can be achieved as long as the pilot keeps the angle of inci- 
dence within an acceptable range. On the other hand if L’(a) is negative 
(i.e. if the lift decreases with increasing angle of incidence) the vibrations 
will build up and may tear the wing apart. 


Exercise 


In designing systems that are capable of vibration it is often desirable to 
introduce friction into the system deliberately, to make the system return 
to equilibrium as quickly as possible after a disturbance. 


If the system is modelled by the equation mX" + kX’ + sX — 0, the prob- 
lem is to find the value of the friction constant k that gives the most rapidly 
decreasing solution for fixed m and s. The rate of decay is indicated by 
the term e~* in the formula for the solution, with the smallest value of a. 
(The larger a, the faster the decay; so we take the smallest value of « as the 
measure of decay, knowing that to be the slowest.) 


Fill in the blanks in the following argument to show that decay is most 
rapid when there is critical damping, i.e. when k? = 4ms. For an exponen- 
tial term e^", we shall call « its rate of decay. 


For small damping, the rate of decay is (i) —  — t since in 
this case 


0 « k? < (ii) 


the rate of decay is less than (iii) = 


For critical damping the rate of decay is 


(iv) = (v) , as k? = 4s, 
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For large damping there are two exponential terms e^*'', e~*' witha, # a2, 
where —a,, —a, are the roots of the auxiliary polynomial equation 


Wi) Hence 
(P + &%)(p + %2) = (vii) 


= (viii) 


and aa, = (ix) 


The rate of damping depends on the smaller of the terms ,, a, and hence 
is less than (x) 


Of the three cases the highest rate of decay is for (xi) 
damping. 


Solution 


G) z, G) 4ms, (i) 1E (iv) a ; 


(v) T (vi) mp? +kp+s=0, 


(vii) p^ + (e, + o)p + oaz, 


k s 

vii) p?+—p+— 

(vii) p ae *x 
@ = 
m 


6) = 


(xi) critical 
9.1.5 Summary of Section 9.1 


In this section we defined the terms 


auxiliary equation (page K128) ub 

or characteristic equation (page K128) FOROR 
small damping (page C16) e 
critical damping (page C16) * 
large damping (page C17) t 

We introduced the notation 
M (page C7) 

and saw that: 

the function M from the set of all constant-coefficient 

linear differential operators to the set of polynomials is an 

isomorphism preserving addition, multiplication and 

scalar multiplication. 

Techniques 

1. Obtain a basis for the solution space of a constant-coefficient second- 

order equation. Wok 

2. Specify the functions that factors of the form (D — à)" and (D? — 

2aD + à? + b?)" contribute to the basis of the kernel of a constant- um d. 


coefficient operator. 
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9.2 EQUATIONS WITH VARIABLE COEFFICIENTS 
9.2.0 Introduction 


If the coefficients in the differential equation are not constant functions, 
then it is normally much harder to solve than the constant-coefficient type 
we considered in Section 1 of this urit. The general theory still applies, 
telling us that if we can find 7 linearly independent solutions of an nth 
order equation then we have a basis for the entire solution space; but these 
solutions are no longer just products of polynomial, trigonometric and 
exponential functions, and indeed they usually cannot be expressed in 
terms of elementary functions at all. In this section we obtain some general 
results that can be used to help overcome this kind of difficulty. These 
general results do not by themselves lead to a complete solution of the 
equation, but they do enable us to recognize a basis of the solution space 
when we see one, and sometimes to complete a basis when we have only 
some of its elements. 


To recognize that a given set of functions, say (yi, ya, ..., Ya}, is a basis 
for the solution space of a given homogeneous differential equation Ly — 0, 
three steps must be taken. We must verify (i) that the given functions are 
in the solution space, (ii) that their number is equal to the dimension of the 
solution space and (iii) that they are linearly independent. These three con- 
ditions are sufficient to ensure that (y,, ..., y,} is a basis, since we know 
that in a vector space of dimension n any n linearly independent vectors 
constitute a basis. The verification of the first of these three conditions is 
simply a matter of checking that the given functions satisfy the equation. 
In this section we will discuss how to verify the second and third. 


We begin with the second condition. Here, it is a question of finding the 
dimension of the solution space. We have already used the fact that this 
dimension is equal to the order of the differential operator L; now we shall 
go into the reason why this is so. 


9.2.1 First-order Equations 


We begin with the simplest case: linear differential equations of order 1. 
We want to show that the solution space has dimension 1. We do this by 
mapping the solution space to another space which we know to have 
dimension 1, namely the space R. The solution space of such an equation is 
a family of functions, whose graphs can be sketched in the manner illus- 
trated below (see Unit M 100 24, Differential Equations I). In the diagram 
the black lines are tangents to the solution curves for y’ = 2y. Typical 
solution curves are shown in red. 
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+ + 
One way of mapping this family of curves to R is to map each curve to its 
intercept on the y-axis. More generally, we could map each curve to its 
intercept on any line parallel to the y-axis; say the line whose equation is 
X= xo With xg some real number. 


In symbols this mapping, which we denote by E, is 
E: y> y(%o) (es) 
where S is the solution set of the first-order linear differential equation. 


It is not difficult to show that E is a linear transformation. (See Exercise 
2 below.) £ is called an initial-condition mapping. 


To complete our argument we want to show that the linear transformation 
E is an isomorphism, i.e. that it is one-one and onto, or equivalently that 
it has an inverse. In other words we want to show that to every real number, 
Say yo, there corresponds a unique function y in the solution set such that 
J(xo) = yo. To do this we need a theorem, which is the subject of the next 
reading passage. 


READ from the beginning of Section 3-4, page K102 to the end of the 
statement of Theorem 3-1 page K103. 


23 


LM 9.2.1 


Notes 


(i) 


Gi) 


(iii) 


Gv) 


(vy) 


line —6, page K102 The equation under discussion is 
aly a) + aolxdyG) =k) (xer) 


which we studied in Unit 4, Differential Equations I. For the present you may 
take it that 4 is the zero function, but the results for other functions A will 
be useful later. 

Equation (3-30), page K102 This was obtained in the preceding section of 
K (pages K96-97). You met it in Unit 4. The expression simplifies con- 
siderably since in our case A = 0. 

line I, page K103 “itis easy to see . . .". In case you do not find it so easy, 
we give the formula for the solution curve (with 4 —0) which passes 
through (xo, yo); it is 


* aox) ) 
= yo exp| — dx 
y —yo o( ee 
(Notice that when x = xo, the integral is zero, and so we get' y = yo. 
Notice also the typical abuse of notation. The xs in aol) dx are dummy 


a(x) 
variables for the purposes of integration and could be replaced by any 


" 


other letter, whereas the x in f is the domain variable in J.) This shows 


E 
that a solution exists; it remains to show that this solution is unique. 

line 9, page K103 “governed by the equation” means here that the particle 
always moves in the direction whose slope is the value of y’(x) given by the 
differential equation, i.e. the slope is (with A = 0) 


__ avo» 
a(x) 


line —10, page K103 “not difficult to show” In case you want to see a 

proof of this, here is one. We follow the suggestions in Exercise 14, page 

K105. 

(a) Lety, and yz be two solutions of the equation: then subtracting using 
the formula (3-30) for the solution, we have 


i69) — y(x) = (ea — c1) exp (- J ails) ax) E 
a(x) 
The exponential part is never zero, so either c; = cı, in which case 
Yi — Ja =0 everywhere in J, or c; # ci, in which case yi—yai*0 
everywhere in J. 
(b) Suppose now that we have an initial-value problem, and suppose that 
we have solutions y, and y; (assumed distinct). Then 


Xi(xo) = yo 
Ya(Xo) 7 Yo 


and so at xo € I, y; — y; is zero. Hence by part (a), y, — y; is zero 
everywhere on J, and the solutions are not distinct. 


Example 


Consider the equation y' = y. 


Then the solution space is S = (x — ce*, c eR} 


We define a mapping E: S —— R, by 


yr— X0) es) 


ie. (x — ce) — ce! = ce, 


This mapping is an isomorphism of the additive Structures. 


(i) 
(ii) 
(iii) 
Gv) 
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What are the images of x ———— &*, x — — 2e* under E? 

What are the inverse images of 1, 2? 

What is the inverse mapping E^! : R — $9 

pis is the solution of the equation Y' = y, with initial condition 
yO) =7? 
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Solution 


@ 


(i) 


(iii) 


Gv) 


Exercises 


By definition E : (x ——— e*) + e! =e 
E: (x — 2e*) + 2e! = 2e 


Since E sends x ~——> e to e, it sends "- to 1. 
e 


1 i 
But 69 €*^U; hence the inverse image of 1 is 


x i— —À e07U and similarly the inverse image of 2 is 
xr 26070. 

Extending (ii) we see that the inverse image of the number r 
is the function x —+ re^? and hence the inverse map- 
ping E^! is given by 

E3ire0—— (x red) — (re R) 

What does the mapping E^! do? Well, it just gives us the 
solution of the equation with a given initial condition at 


1. Hence if we want a solution with y(1) = 7, then all we need 
is E71 (7) = x ——25 Te), 


1. S denotes the solution set of the differential equation 


yQ)-2xyx) (eR) 


and E is the function from S to R defined by 


Ey ——5 yQ). 


(i) Does the function x ——— exp (x?) belong to S? 
(i) Write down the image of this function under E. 
(iii) Write down the inverse image of the number exp (4) under E. 


(v) If 


yo denotes an arbitrary real number, write down its inverse 


image under E. 
(v) Write down the inverse of the function E. 
2. Verify that E, as defined in the text on page C23, is a linear transfor- 


mation. 


Solutions 


1. 


() Yes. 

(i) The image of y under E is defined to be (2); so if 
y(x) = exp (x?) the image of this y is exp(4). 

(iii) Since E maps the function x ————» exp (x?) to exp (4), 
the inverse image of exp (4) is this function. 

(iv) The inverse image of the number yp is the solution of 
the differential equation such that y(2) = yg. Since the 
general solution is y(x) = c exp (x?), the required value 
of c is determined by yo = c exp (27), i.e. c = yg e^ * and 
so the required solution is given by 


V(x) = yo exp (x? —4) — (xeR) 
The inverse image of yo is thus the function 
x= yoexp(x?-4) (xe R). 
(V) E`: yoi Qo yo exp (x? —4) — (Yo ER). 


To show that E is a linear transformation we must show that 
its domain is a vector space and that 


E(ay + Bz) = aE(y) + BE(z). 


The domain, S, is a vector space since it is the solution space 
of a homogeneous linear problem Ly = 0. 
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The transformation is linear since 


E(ay + Bz) = (ay + Bz)(xo) (definition of E) 
= ay(Xo) + fz(xo) (definitions of addi- 
tion of functions 
and scalar multi- 
plication) 
= aE(y) + BE(z) (definition of E) 


9.2.2 Equations of Higher Order 


To illustrate how the ideas explained in the last sub-section can be extended 
to higher-order equations, let us consider a second-order equation, say 


=x) (xen a) 
We hope to show that its solution space S is two-dimensional, i.e. that it 
is isomorphic to R?. One way of setting up the mapping from S to R? is to 
rewrite the equation as a pair of simultaneous first-order equations, in 
terms of y and a new function z = y’. 

z(x) = xy) 

h xe 2) 

y&-z CED : 

The solution set of this pair of equations is a set of pairs of functions of the 


form (y, z), and it is equivalent to the solution set S of Equation (1) in the 
sense that y e S if and only if the pair (y, y") is a solution of Equation (2). 


The natural generalization of the initial-condition mapping E, which we 
looked at in the previous sub-section, to the System of simultaneous 
equations (2) is the mapping 


O, z) ——9 (yo), z(xo)) — ((v, z) € solution set of Equation (2)) 
where xo e I. This mapping is equivalent to 


0, y) ——9 (y(%o), yx) Wes). 


This idea leads to a natural way of mapping the solution space S of our 
original Equation (1) to R?; as before we shall denote this initial-condition 
mapping by E: 


E: y — (yo), Y(xQ) Qes) 


Once again, it is not hard to show that E is a linear transformation. To 
complete our proof that S is two-dimensional, we want to show that Eis 
an isomorphism, i.e. that it has an inverse. This is done by proving that to 
every ordered pair of real numbers (yo, Yo) there corresponds a unique 
function y in the solution set such that (xo) = yo and y'(x,) = y4. The 
theorem which proves this is stated (but not proved) in the next reading 
passage. 


READ from line —6, page K103 to line 13, page K104. 


Notes 


() line —1, page K103 The proof referred to will be treated in Unit 33, 
Existence and Uniqueness Theorem for Differential Equations. 

(ii) Equation (3-32), page K104 The notation is a little different from ours: 
the book uses y, for our y4, ya for our yb, etc. 


The main result of the last two sub-sections can be summed up in the 
following theorem: 
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Theorem 


If L is a normal linear differential operator of order n, with domain C'(I), 
S is the solution set of Ly = 0, and xo is any point in J, then the mapping 
E from S to R" defined by 


Ei y ——95 (69). Y (9), ..., y 7 (xg) 
is an isomorphism. 


Example 
Consider the differential equation 
»'(3) * y) 20 
and initial point xg — 0. The equation has solution space S with basis 
{x ——> l, x +e}. Let E: S ——+ R? be the initial-condition 
mapping 
E: y — Wo) Y G9) Wes). 
(i) What is 
(a) E(x —— 1) 
(b) E(x——— e™), 
() E^Q,0) 
(d) E^!(0,1)? , 
(i) What is E !(c, c2)? 


Solution 
@ (a). If y:x — — 1, then y’:x ——— 0 
E(y) = Wo), y'(xo)) = Q0), y) = (1, 0) 
(b) If yix — 5 e™*, then y: xm —e^* 
E(y) = (vO), »'(0)) — (1, — 0. 
(c) From (a) we have E^!(1, 0) = x ——+ 1 


(d) Wecan use the results in (a) and (b). The pairs (1,0) and 
(1, — I) are a basis for R? and 


©, 1) — (1,0) - (1, —1). 


Using the fact that E~', an isomorphism, is a linear 
transformation (see Theorem 1.1, page N28) we have 


E~'(0, 1) = £11, 0) — E“(1, —1) 
= (x —> D)-(x——e7) 
-—Xx-———l1—e* 
(i) We can use the same technique as in (i) (d). 
(e, 62) = (e + €2)(1, 0) — e2(1, — 1) 
so that E^!(&, c2) = x ——9 ¢, + e — ce. 
The function E~'(c;, c3) is just the solution of the original 
equation y'(x) + y'(x) =0 which satisfies the initial condi- 
tions y(0) = c, y'(0) = c3. 


Exercises 


l. For the differential equation 
V@)+ye)=0 (ER) 
with xo =0, write down E(x.——— cos x), E(x ——— sin x), 


E~'(1, 0), E74(0, 1) and E~(c,, c;). Note that the solution space has 
basis {cos, sin}. 
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2. For the differential equation 
y'G) — y099 = 0 
with xg =0, write down E(x ——> e*), E(x ——5 e?) EM, 0), 
E~*(0, 1), E^! (c, cz). Note that the solution space has basis 
(x —— e, x c e. 
Exercise 2, page K105. 
4. Exercise 15, page K105. 


Ld 


Solutions 


1. E(x:.——— cos x) = (cos 0, cos’ 0) 
= (1, 0), since cos’ = —sin 
E(x ———— sin x) = (sin 0, sin’ 0) 
= (0, 1) since sin’ = cos 
E^ (1, 0) = x ——— cos x 
E^ X0, 1) = x —— sin x 
E7*(¢,, c) = x ——— c, cos x + c; sin x 


We can represent this solution by a diagram as in the T.V. 
programme 


R? 
(1,0) 


(cosO, cos‘0) 


(0,1) 
| »(sinO, sin’O) 


y |—-—E.—s (yoyo) 


2. E(x — —» e”) = (exp 0, exp’ 0) 
=(1,)) 
E(x ——— e™*) = (1, -1) 
E“(1, 0) = E Q1, 1) +40, —1) 
= x> }(e* + e7”) 
E“(0, 1) = EQU, 1) — 4, - 1) 
=x Ke - e) 
E^ (6, 62) = x ——95 ee" + e") + dex" — e~”) 
3. If we want the equation to be normal, we must choose an 
interval in which sin x does not vanish and which includes 
the initial value = Such an interval is (0, z) = (x:0 < x <x}. 


In this interval we can apply standard techniques (see Unit 4, 
Differential Equations I) to obtain the general solution in the 
form 


y(x) = efsin x. 


Using the initial condition, we get 
3 = 
€ =2= c2 
4 
c= J2 


Thus, the solution to the initial value problem is 


yx) = of 2jsin x (x e (0, x). 


4. A solution to the exercise is given on page K733. Alternatively 
we could use the equation in Exercise 3 above. It is not normal 
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on the interval (0, zt) = (x: 0 <x < z}, and we cannot find a 
solution for which y(0) 4 0, because if we put x — 0 in the 
equation, we get 


(sin 0)y(0) + (cos 0)y(0) = 0 
ie. x0) = 0. 
So y(0) cannot be specified arbitrarily. 


9.2.3 Dimension of the Solution Space 


The theorem stated at the end of the last sub-section gives the dimension 
of the solution space of a normal nth-order linear differential equation; 
for this solution space is isomorphic to R", which has dimension n, and 
isomorphic spaces have equal dimensions. This result is consolidated in the 
next reading passage. 


READ Section 3-5 of K, omitting the proof of Theorem 3-3 and the follow- 
ing paragraph. 


Notes 


(i) line —7, page K106 The C(I) in the statement of the theorem would have 
been better as C"(I) see page K92. It is not wrong, however, because C"(7) 
is a subspace of C(Z). Over the next few pages the authors skip about a bit: 
it is always best to use the most appropriate subspace; “ most appropriate” 
is determined by the highest derivative required in the context. 

(ii) line —6, page K106 You need not read this proof, since the theorem fol- 
lows immediately from the existence of the isomorphism E, as noted above. 

(iii) Example 1, page K108 In this example, the standard basis in Æ? is 
mapped into S and yields a basis in S. The functions forming this basis, 
denoted by cosh and sinh, are called the Ayperbolic cosine and the hyper- 
bolic sine.* You have already met them in Exercise 2 of the preceding sub- 
section, but they were not introduced properly. They have properties 
analogous to those of the cosine and sine functions; this analogy arises 
from the similarity between the formulas defining cosh and sinh and for- 
mulas for cos and sin which can be derived from Euler’s formula (see 
Unit M100 29, Complex Numbers II). 


1 
cosh x = I (e*+e7*), cosx— 2 (e'* + e7?) 


1 _! X =x " ui ft La tx 
sinh x = 5 (e —e7?), sinx = 5; (e e7"). 


As an example of this analogy, we saw in Exercises 1 and 2 of the preceding 
sub-section that for the equation y” + y = 0, we have E~'(1, 0) = cos and 
E-(0,1) =sin; but, for y"—y —0, we have E-'(1,0) — cosh and 
E-*(0, 1) ^ sinh. 


4 


4 


Xi—cosh x 


x— sinh x 


* The functions cosh and sinh are usually read as “cosh” and “shine”, respectively. 
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(v) Example 2, page K109 This time we use ((1, 1), (1, —1)) as our basis in 
4?*, and the corresponding basis in S is {e*, e-*}. 

(v) Corollary 3-1, page K109 In terms of the isomorphism E, this corollary 
may be stated as follows: if B is a basis for g", then E-'(B) is a basis for S. 


Example 


(i) Show that the functions x-——+ cosh x, x ———»^ sinh x and 
X > cosh (a + x), where a is a constant, are all solutions of 


(D? - y 20. 


(ii) By considering E: S +. R? for initial point x9 — 0, find a linear 
relation between them. 


Solution 


0) 


(ii) 


Exercises 


l. (i) Show that x—— sin (a + x), 


From the formulas 
cosh x = 1(e* + e^?) 
sinh x = à(e* — e~*), 
we find that cosh'(x) = sinh x and sinh'(x) = cosh x. Hence 
cosh'(x) = sinh'(x) = cosh x 
cosh'(x) — cosh x = 0 
(D? — 1) cosh = 0. 
Similarly, (D? — 1) sinh =0 


Since ais a constant, the function x —-—» cosh (a + x) also 
satisfies the differential equation. 


E: S ——+ R? 
is specified byE: y —> (y(0), y'(0)) 
Hence 
E(x —> cosh x) = (1, 0) 
E(x — sinh x) = (0, 1) 
E(x — cosh (a + x)) = (cosh a, sinh a) 
Since E is an isomorphism, any linear relationship between 
E(x —> cosh x), 
E(x ———— sinh x), 
E(x —> cosh (a + x)) 


corresponds to one between the functions cosh x, sinh x, 
cosh (a + x). 


Now, 
E(x +> cosh (a + x)) 
$ = cosh aE(x > cosh x) 
+ sinh aE(x ——> sinh x). 


Hence, cosh (a 4- x) = cosh a cosh x + sinh a sinh 


X 0———— sin x, X —— —5 cos x 


satisfy (D? + 1)y = 0. 
(ii) By considering the initial-condition mapping E: S ——— R? 
with xo = 0, find a linear relation between them. 


2. (i) Show that the functions x:—— I, 


xe cos 2x, and 


X > cos? x are solutions of (D? + 4D)y =0. 


(ii) By considering the isomorphism E: S ——+ R? 
find any linear relations between them. 
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with xo — 0, 


Solutions 


1. (i) Differentiation shows that they satisfy the equation. 


(i) 


(ii) 


Their images under E are 

E(x +——> sin x) = (sin 0, sin’ 0) = (0, 1) 

E(x ———» cos x) = (cos 0, cos’ 0) = (1, 0) 

E(x +——> sin (a + x)) = (sin a, sin’ a) 

= (sin a, cos a) 

In R? we have the linear relation 

(sin a, cos a) = (sin a)(1, 0) + (cos aY(0, 1). 
Applying E^! to both sides and using the fact that it is 
in an isomorphism we find 


E^! (sin a, cos a) = sin a E^! (1, 0) 
+ cosa E^ (0, 1) 


sin (a + x) = sina cos x + cos a sin x. 


Differentiation shows that they are solutions of the 
equation. Their images under E are 

E(x ——— — 1) = (1,0, 0) 

E (x +—— cos 2x) = (cos 0, —2 sin 0, — 4 cos 0) 

= (1,0, -4) 
E (x ———5 cos ?x) = (cos? 0, —2 cos 0 sin 0, 
—2 cos?0 + 2 sin? 0) = (1, 0, —2). 

To find the linear relation between these vectors we can 
use inspection or the method of Unit 3, Hermite Normal 
Form, (pages N72 and N73). Using the latter method we 


consider the matrix whose columns are the coordinates 
of the 3 vectors, 


H^ dq. ou 
0 0 ol, 
[0 -4 -2 


and by row operations reduce it to its Hermite normal 
form, which is 


[10g 
01 4 
lo 0 0 


The last column is half of the first plus half of the second 
and so the same relation holds for the columns of the 
first matrix 


(1, 0, —2) = 4(1, 0, 0) + 3(1, 0, —4). 
Applying E^! to both sides of the equation we obtain 


cos? x = 1 (cos 2x + 1) 


a relation which is useful, for example, if you want to 
integrate cos? x. 
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9.2.4 The Wronskian 


For the general solution of a linear homogeneous equation it is not enough 
to have as many solutions as there are dimensions in the solution space; 
we must also make sure that these solutions are linearly independent, so 
that they form a basis. For example, in the last exercise of the preceding 
section, three solutions of the equation (D? + 4D)y = 0 are the functions 
X ——— l, x> cos2x and x+——>+ cos? x, but, since they are 
connected by the linear relation 1 + cos 2x — 2 cos? x = 0, they are not 
linearly independent, and so cannot be taken as a basis. In this sub- 
section we generalize the idea of that exercise, using it to set up a general 
theoretical criterion for the linear independence of functions. 


The criterion we shall set up is based on the determinant function, which 
you studied in Unit 5, Determinants and Eigenvalues. Remember that the 
determinant of any n vectors a,,..., a, in R” is a scalar D(a,, ..., 8,), 
whose main property is that it vanishes if and only if the vectors are linearly 
dependent. We can also write the coordinates of a}, ... , a, (with respect to 
the standard basis in R") as the columns of an n x n matrix, say A, and then 
D(a,, ..., a,) is called the determinant of that matrix, written det A. For 
example, in R?, if a, — (p, q) and a; — (r, s) then we write 


pr 


D(a,, a2) = det [ E 5 s 


qs 

= ps — rq. 
Let us apply this idea to the three solutions of (D? + 4D)y = 0 just con- 
sidered. We saw in the exercise that the images of the three solutions 
x— l, x > cos 2x, x —— —» cos? x under the isomorphism E 
are (1, 0, 0), (1, 0, —4) and (1, 0, —2). The corresponding determinant is 


l 1 1 Io. å 
det|O 0 Ol,ielo 0 O0 
0 -4 -2 0 -4 -2 


and it is equal to zero because one of its rows contains only zero entries. 
The vanishing of this determinant tells us that the three image vectors 
(1, 0, 0), (1, 0, — 4) and (1, 0, —2) are linearly dependent. Thus, since E is 
an isomorphismit follows that their inverse images, the functions x —— I, 
X —— —À5 cos 2x and x ++ cos? x are also linearly dependent, and do 
not constitute a basis for the solution space. Thus, if all we want to know is 
whether the functions constitute a basis, we can avoid the explicit calcu- 
lation of the linear relations as we did in the exercise; the determinant tells 
us at once that some linear relation exists, and that is enough to rule out, 
the given functions as a possible basis. On the other hand, if we had found 
the determinant to be different from zero, we would have known that its 
columns were linearly independent, and hence the original functions 
(which were known to span the space) would have been linearly indepen- 
dent and hence would have been suitable as a basis of the solution space. 


Y 


Generalizing these considerations we have the following result: 


Theorem 1 


A set of n solutions y,, ..., y, of an nth-order normal homogeneous linear 

differential equation is linearly independent if and only if the determinant 
DEY), ..., E(y9) 

is not zero. 

The proof is a consequence of the fact that E is an isomorphism: the func- 

tions yi, ..., y, are linearly independent if and only if their images under E 


are linearly independent, and these are linearly independent if and only if 
the determinant of this set of image vectors is not zero. 
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There is another way of stating this theorem, which makes use of the fact 
that the mapping 


E: y — (y(X0); Y Gio), «+ 9 7G) 


depends on the number xo. This has the consequence that the determinant 
mentioned in the theorem depends on x, ; written out as the determinant 
of a matrix it looks like this: 


Ji (xo) "SS Xen) 
DEO), -s EO) = det | 9) — o) 


EMR) oo IE) 
This means that we can define a function (with domain 7 and codomain R): 
W : xo —> (the above determinant) (xo EI). 


It is called the Wronskian function associated with the set of functions 
Pis +++ Yn: 


In terms of the Wronskian function we can re-state the above theorem in 
function language. 


Theorem 2 


A set of n solutions of an nth-order normal homogeneous linear differential 
equation on the interval I is linearly independent if and only if its 
Wronkskian function maps all the points in 7 to numbers other than 
zero; and it is linearly dependent if and only if its Wronskian function is 
the zero function. 


Section K3-6 gives some examples of the use of the Wronskian, beginning 
with cases where y,,..., y, are not necessarily solutions of a differential 
equation. In this case E is no longer an isomorphism, and so less can be 
proved; if the Wronskian function is ever non-zero, then the functions are 
linearly independent, but if it is the zero function, we can deduce nothing. 
We have not included any of this in the official reading passage, since we 
are mainly interested in the case where y,, ..., y, do satisfy a differential 
equation. 


READ (a) page K111 down to the equation “ W[x, 2x] +: — 0"; (b) the 
statements only of Theorems 3-5 and 3-6 on page K113; (c) the whole of 
page K114. 


Notes 


(i) Equation (3-41), page K111 These vectors are just E( y,), ..., E( 4). 
(i) line 7, page K114 “linearly independent in C(7)” that is, the functions 
sin? x and 1/sin? x, both of which belong to (J), are linearly independent. 


Exercises 


1. Examples 4 and 5, page K114 refer to an interval Jin which tan x and 
cot x are both defined. Give an example of such an interval. 


2. Use Theorem 2 to test the following sets of functions for linear inde- 
pendence, and write down a basis for the solution space if there is 
one among the functions given. In each case, take 7 to be R. 


(i) ee (consider (D? — 1)y = 0) 
(ii) sin x, cos x (consider (D? *0Dp20) 
(i) 1, x, x? (consider D*y = 0) 


(i) 1-xx-x*,x!-1 (consider D?y = 0) 
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Solutions 


1. Any interval that includes no integer multiple of 5 will do, 


for example, Q 5) (vu not fo, jJ 


: sinx ,.. 3 
Since tan x = ——, it is defined unless cos x = 0, i.e. unless 
cos x 
A X ln cosx .., 
x is an odd multiple of y Since cot x =, it is defined 
sin x 


$ * ui " 
unless x is an even multiple of 5. Thus the interval must not 


include either even or odd multiples of 5. 
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() Both functions satisfy the given differential equation, and 


le" e" 


War. Sa 


= 7250; 


so they are linearly independent. Since n —2 here, 
{e*, e~*} is our basis. 
(ii) Both functions satisfy the differential equation, and 


sin x cos x] 
cosx  —sinx 


= —sin? x — cos? x 
=-1#0; 


so they are linearly independent, and since n =2, a 
basis is {sin x, cos x}. 
(iii) All three functions satisfy the equation, and 


1 x x 
0 1 2x 
002 


W= =2 


since W is triangular. 


Hence, the functions are linearly independent and since 
the equation is third-order, a basis is (1, x, x°}. (In a 
similar way we can show that {1, x, x?, ..., x?) is 
linearly independent, for any positive integer 7.) 

(iv) All three functions satisfy the equation, and 


2 


1-x x-x? x*-1 


W=| -1 1-2x 2x 

0 -2 2 

1-x x-1 x!-1 
-|-l 1 2x| (col. 2 + col. 3) 

0 0 2 


= 0, since first 2 columns are proportional. 


The given functions are linearly related, and do not form 
a basis. 
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9.2.5 Summary of Section 9.2 


In this section we defined the terms 


initial-condition mapping (page C23) RN 
initial-value problem (page K103) OX 
hyperbolic sine—sinh (page C29) EE 
hyperbolic cosine—cosh (page C29) abd 
uniqueness problem (page K103) kea 
Wronskian function (page C33) bes 
We introduced the notation 
E (page C23) 
Theorems 
1. (3.1, page K103) 
Every initial-value problem involving a normal first-order linear differential E Or E 
equation has precisely one solution. 
2. (3.2, page K104) (The existence and uniqueness theorem for linear dif- xc 
ferential equations) 
Let 
aya) 52, asy m AG) A 


be a normal nth order linear differential equation defined on an interval 
I, and let x, be any point in Z. Then if y, ..., y, ., are arbitrary real 
numbers, there exists one and only one solution y(x) of Equation (A) with 
the property that 


J(Xo) = Yo, Y (xo) = Ju «+s y YG) nS 
3. (page C27) oro 
If Lis a normal linear differential operator of order n, with domain cn, S 


is the solution set of Ly = 0, and x, is any point in 7, then the mapping E 
from S to R” defined by 


E: y —> (yo), Y G9), ..., y 79) 
is an isomorphism. 


4. (3-3, page K106) . 
The solution space of a normal nth order linear homogeneous differential * oko 
equation is z-dimensional. 


5. (Theorem 1, page C32) 

The solution y,,..., y, of an nth-order normal linear homogeneous *ox 
differential equation is linearly independent if and only if the determinant 

D(E(y4), ..., E(y,)) is not zero. 


6. (Theorem 2, page C33) . 

A set of n solutions of an nth-order homogeneous linear differential *ox 
equation on the interval I is linearly independent if and only if its 
Wronskian function maps all the points in J to numbers other than zero; 

and it is linearly dependent if and only if its Wronskian function is the 

zero function. 


Techniques 

1. Use the initial condition mapping E and its inverse E^! to solve initial trx 
value problems. 

2. Use E to determine linear relations between functions in the solution * 
space of a linear differential equation. 

3. Use the Wronskian function to determine if a set of solutions of a * 


linear differential equation forms a basis. 
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9.5 SUMMARY OF THE UNIT 


The general aims of this unit are to show how to obtain the solution space 
for a constant-coefficient linear differential equation, given that the opera- 
tor is factorized into first and second-order factors, and how to recognise a 


basis for the solution space of any homogeneous linear differential equa- 
tion. 


We began by showing that factorization of constant-coefficient differential 
operators is the same as that for polynomials and then showed how the 
solutions of a second order equation depended on whether the roots of the 
corresponding quadratic are real and distinct, real and the same or 
complex conjugates. 


If the equation is (D? + a, D + ag)y = 0 and 
(D? + aD + a9) =(D — a)(D — f) 
then for 


(i) «s f, both real, the general solution is c,e** + c; e°" 
(i) a = B, real, the general solution is (c, + c; x)e** and 
Gii) «,f complex, =a + ib, b #0, the general solution is 
e™(c, cos bx + c; sin bx). 


Then we showed how these results could be extended to an operator of 
any degree, as long as we could factorize it. Typically the factor (D — a)" 
would contribute the functions e**, ..., x' "!e"* to the solution space and 
the one of the form (D? — 2aD + a? + b?y would contribute e°% cos bx, 
«++, x77 1e™ cos bx, e™ sin bx, ..., x" 4e™ sin bx. 


The first section ended with an example of modelling a physical situa- 
tion, an aircraft wing, by differential equations, and discussed the relation 
between the coefficients of the equation and stability or damping. 


In the second section we discussed general homogeneous linear differential 
equations and stated that for an nth order equation there exists a unique 
solution for each initial condition (xo), ..., "^ (xg). In this way we are 
able to set up the initial-condition mapping E from the solution space to R" 
which is an isomorphism. This tells us that the solution space is n-dimen- 
sional and that if we have n functions in it which we wish to test for linear 
independence then we only have to test their images in A". Hence the 
functions y,, ..., Ya are independent if and only if det (E(y;), ..., E(y,)) is 
not equal to zero. If they are linearly dependent then we can use Hermite 
normal form to discover any linear relations between them by finding the 
relations between their image vectors. 


Definitions 
auxiliary equation : (page K128) 

or characteristic equation (page K128) 
small damping (page C16) 
critical damping (page C16) 
large damping (page C17) 
initial-condition mapping (page C23) 
initial-value problem (page K103) 
hyperbolic sine—sinh (page C29) 
hyperbolic cosine—cosh (page C29) 
uniqueness problem (page K103) 
Wronskian function (page C33) 
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* xoxo Xo 


* 
+ +++ 
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Theorems 


1. (3.1, page K103) . 
Every initial-value problem involving a normal first-order linear differen- 
tial equation has precisely one solution. 


2. (3.2, page K104) 
(The existence and uniqueness theorem for linear differential equations) 
Let 


a9) TE + aglady = HO) (A) 


be a normal nth-order linear differential equation defined on an interval 
I, and let x, be any point in J. Then if yo, ..., y,-1 are arbitrary real 
numbers, there exists one and only one solution y(x) of Equation (A) with 
the property that 


X(xo) = Yo; Yo) = Y» oe ye) = PSU 
3. (page C27) 
If L is à normal linear differential operator of order n, with domain 


C'(I), S is the solution set of Ly = 0, and x, is any number in J, then the 
mapping E from S to R" defined by 


E: y — 95 (xo) y G9), -> Yo) 
is an isomorphism. 


4. (3-3, page K106) 
The solution space of a normal nth order linear homogeneous differential 
equation is n-dimensional. 


5. (Theorem 1, page C32) 

The solution y,,..., y, of an nth-order normal linear homogeneous 
differential equation is linearly independent if and only if the deterniinant 
D(E(y4), ..., E(y,)) is not zero. 


6. (Theorem 2, page C33) 

A set of n solutions of an nth-order homogeneous linear differential 
equation on the interval 7 is linearly dependent if and only if its Wronskian 
function maps all the points in 7 to numbers other than zero; and it is 


linearly dependent if and only if its Wronskian function is the zero 
function. ` 


Techniques 


l. Obtain a basis for the solution space of a constant-coefficient second 
order equation. 

2. Specify the functions that factors of the form (D — a)" and (D? — 2a D 
+ a? +b”) contribute to the basis of the kernel of a constant-coefficient 
operator. 

3. Use the initial-condition mapping E and its inverse E^! to solve 
initial-value problems. 

4. Use E to determine linear relations between functions in the solution 
space of a linear differential equation. 

5. Use the Wronskian function to determine if a set of solutions of a 
linear differential equation forms a basis. 


Notation 
M (page C7) 
E (page C23) 
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9.4 SELF-ASSESSMENT 
Self-assessment Test 


This Self-assessment Test is designed to help you test quickly your under- 
standing of the unit. It can also be used, together with the summary of the 
unit for revision. The answers to these questions will be found on the 


next non-facing page. We suggest you complete the whole test before 
looking at the answers. 


1 Which of the general solutions 


(a) Ae7?*4 Be™ 

(b A+ Be^?* 

(c) Axe^?* + Be7?* 

(d) 4e^?* cos 2x + Be^?* sin 2x 


corresponds to the following differential operators: 


(i) D?+4D+8 
(i) D?+4D+4+4 
(iit) D?42D 
(iv) D?-4 
2 Write down a basis for the solution space for the following differential 
equations: 


(i) (D—-1)(D —4)(D + 5)y =0 
(i) (D+ 1) Dy =0 
Gi) (D? + I/D — Dy =0 
Gv) (D' 4D 5)D!y-0 

() (D? -4D 45! (D—1y =0 
(vi) (D? + 1D! 499y—0 


3 Which of these functions occur as solutions to a homogeneous con- 
stant-coefficient second-order linear differential equation? 


(i) x= e 
(i) x:—— cos x 
(iii) x ——> xe?* 
(i) x—— — xsinx 
() x— xe 
(vi) x ——— cos 2x 
(vii) x —— cos? x 


4 Find the simplest constant-coefficient linear differential operator which 
annihilates the following functions 


(i) x———9x-ccosx —3sinx 
(i) x e +3 + 4xe* 
(ii) x ——+ (x + Ie — e7”) 
5 What is the dimension of the solution space of 
() D?y=0 
(i) (D? + 3D)y =0 
(i) (D+2)y=07 
6 Associated with the solution space S = (y: (D — 1)!y « 0) and the 
initial point 0 we have the initial-condition mapping E: S ——> R?, 
What is 
G E@r—e*) 
(i) E(x——— xe”) 
(Gu) E^!(1,0)? 
7 Given that cos 3x, cos x, cos? x are solutions of (D? + 9)(D? + l)y = 0, 
find any linear relations between them. 
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Solutions to Self-assessment Test 


1. We check which general solution corresponds to which differential 
operator either by seeing if the operator annihilates it or alternatively by 
factorizing the operator, constructing its general solution and then looking 
for it in the list. 


() D?+4D+8=D? —2aD +a? +b? where a= —2,b — 2. Hence 
the basis of its solution space is (e ?* cos 2x, e ?* sin 2x} and the 
general solution is (d). 

(ii) D?+4D+4=(D+2)*. The auxiliary polynomial has repeated 
root —2 and hence a basis for its solution space is {e~?*, xe~?*}. 
This corresponds to general solution (c). 

(iii) D? +2D = D(D + 2). The auxiliary polynomial has roots 0, —2. 
The solution space has basis (1, e~?*} and the general solution is (b). 

(iv) D*—4-(D-—2Y(D-2) The auxiliary polynomial has roots 2, 
—2 and hence the solution space has basis (e?*, e~?*} corresponding 
to general solution (a). 


2. Since all the operators are already factorized as much as possible over 
the reals we can just write down the bases 


à) (eee 
(i) {e7*, xe^*, x7e7*, 1, x) 
(iii) (cos x, sin x, e?) 
(iv) {e?* cos x, e?* sin x, 1, x} 

(v) (e* cos x, xe?* cos x, x?e?* cos x, e?” sin x, xe?" sin x, x?e?" sin x, e") 
(vi) {cos x, sin x, cos 3x, x cos 3x, sin 3x, x sin 3x). 


3. (i) Yes: it is a solution of (D — 1)Dy = 0 
(i) Yes: it is a solution of (D? + Dy =0 
(ii) Yes: it is a solution of (D — 27y = 0 
(iv) No: as the lowest order homogeneous constant-coefficient 
equation it satisfies is (D? + 1)?y = 0, fourth order. 
(v) No: as the lowest order homogeneous constant-coefficient 
equation it satisfies is (D — 1)?y = 0, third order. 
(vi) Yes: it is a solution of (D? + 4)y — 0. 
(vii) No: cos? x = 4(1 + cos 2x) and hence the lowest order homo- 
geneous constant-coefficient equation it satisfies is D(D? + 4)y 
= 0, third order. 


4. (i) x is annihilated by D?,cos x and sinx by (D? +1). Hence 
D?(D? + 1) is the simplest annihilator of x + cos x — 3 sin x. 

(ii) 3isannihilated by D, e" and xe* by (D — 1)?. Hence D(D — 1)? 
is the simplest annihilator of e* + 3 + 4x e. 

(iii) e", xe* are annihilated by (D — 1)? and e^", xe^* are annihilated 
by (D41Y. Hence (D- D(D - 1? 2(D? 1? is the 
simplest annihilator of (x + 1)(e* — e7*). 

5. The dimension of the solution space is just the order of the equation, so 
we have 


(i) 2, (i) 2, (i) 1. 
6. The initial-condition mapping for initial point 0 is just E: S + R 
y —> (0), YM) (eS). 
Hence 
(G) E(x —— e) = (e°, e?) = (1, 1) 
Gi) E(x — — xe") = (0e, (0 + 1e?) = (0, 1) 
and 
(i) E-(1,0)2 E-(1, 1) — (0, 1) 
= E“'(1, 1) — E7*(0, 1) 
= et — xe* =(1 — x)e*. 
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7. Given that cos 3x, cos x, cos? x are solutions of the fourth-order 
equation we consider the initial-condition mapping at the point 0: 
E: S — R* 
3 — — 90), YO, YO, y") (yes), 

where S is the solution space of the equation (D? + 9(D? + 1)y =0. 
We find E(x ——> cos 3x) = (1, 0, —9, 0) 

E(x —— 5 cos x) = (1,0, —1,0) 

E(x —— — cos? x) = (1,0, —3, 0) 


as (cos*)'x = —3 cos? x sin x (cos?)'(0) = 0 
(cos?)"x = 6 cos x sin? x — 3 cos? x (cos?)'(0) = —3 
(cos?)"x = —6 sin? x + 12 cos? x sin x + 9 cos? x sin x (cos?)"(0) — 0 


Since E is an isomorphism any linear relation between E(cos 3x), E(cos x) 
and E(cos? x) corresponds to a linear relation between cos 3x, cos x and 
Cos? x. It is easy to see that we only need consider relations between the 


EL 


We reduce to Hermite normal form the array 


1 1 1 
-1 -3 -9 


1 0 -3 
0 1 4 
Hence 

1 1 1 
0 0 0 

am tes eg 
0 0 0 

E (x — — cos 3x) = 4E (x > Cos? x) — 3E (x +> cos x) 


and since E is an isomorphism we have 


cos 3x = 4 cos? x — 3 cos x. 
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LINEAR MATHEMATICS 
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Vector Spaces 

Linear Transformations 

Hermite Normal Form 

Differential Equations I 

Determinants and Eigenvalues 

NO TEXT 

Introduction to Numerical Mathematics: Recurrence Relations 
Numerical Solution of Simultaneous Algebraic Equations 
Differential Equations II: Homogeneous Equations 
Jordan Normal Form 

Differential Equations III: Nonhomogeneous Equations 
Linear Functionals and Duality 

Systems of Differential Equations 

Bilinear and Quadratic Forms 

Affine Geometry and Convex Cones 

Euclidean Spaces I: Inner Products 

NO TEXT 

Linear Programming 

Least-squares Approximation 

Euclidean Spaces II: Convergence and Bases 

Numerical Solution of Differential Equations 

Fourier Series 

The Wave Equation 

Orthogonal and Symmetric Transformations 
Boundary-value Problems 

NO TEXT 

Chebyshev Approximation 

Theory of Games 

Laplace Transforms 

Numerical Solution of Eigenvalue Problems 

Fourier Transforms 

The Heat Conduction Equation 

Existence and Uniqueness Theorem for Differential Equations 
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